2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method          下载免费PDF全文

Yuan Yuan  Jianke Qiang  Jingtian Tang  Zhengyong Ren  Xiao Xiao 《Geophysical Prospecting》2016,64(3):767-779

To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion.  相似文献   

Seismic design and testing of buckling‐restrained braces with a thin profile          下载免费PDF全文

Pao‐Chun Lin  Keh‐Chyuan Tsai  Chieh‐An Chang  Yu‐Yun Hsiao  An‐Chien Wu 《地震工程与结构动力学》2016,45(3):339-358

A thin‐profile buckling‐restrained brace (thin‐BRB) consists of a rectangular steel casing and a flat steel core that is parallel to a gusset plate. A thin configuration reduces the width of the restraining member and thus saves usable space in buildings. However, deformable debonding layers, which cover the steel core plate in order to mitigate the difference between the peak tensile and compressive axial forces, provide a space for the steel core to form high mode buckling waves when the thin‐BRB is under compression. The wave crests squeeze the debonding layers and produce outward forces on the inner surface of the restraining member. If the restraining member is too weak in sustaining the outward forces, local bulging failure occurs and the thin‐BRB loses its compression capacity immediately. In order to investigate local bulging behavior, a total of 22 thin‐BRB specimens with a ratio of steel core plate to restraining steel tube depth ranging from 0.3 to 0.7 and axial yield force capacities ranging from 421 kN to 3036 kN were tested by applying either cyclically increasing, decreasing, or constant axial strains. The restraining steel tube widths of all the specimens were smaller than 200 mm and were infilled with mortar with a compressive strength of 97 MPa or 55 MPa. Thirteen of the 22 thin‐BRB specimens' restraining members bulged out when the compressive core strains exceeded 0.03. A seismic design method of the thin‐BRB in preventing local bulging failure is proposed in this study. Test and finite element model (FEM) analysis results suggest that the outward forces can be estimated according to the BRB compressive strength, steel core high mode buckling wavelength, and the debonding layer thickness. In addition, the capacity of the restraining member in resisting the outward forces can be estimated by using the upper bound theory in plastic analysis. Both the FEM analysis and test results indicate that the proposed method is effective in predicting the possibility of local bulging failure. Test results indicate that the proposed design method is conservative for thin‐BRB specimens with a large steel core plate to restraining steel tube depth ratio. This paper concludes with design recommendations for thin‐BRBs for severe seismic services. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

基于二次场电导率分块连续变化的三维可控源电磁有限元数值模拟   总被引：5，自引：5，他引：0       下载免费PDF全文

李勇  吴小平  林品荣 《地球物理学报》2015,58(3):1072-1087

从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态.  相似文献   

台湾碰撞带现今地壳变形场特征及其动力学成因的有限单元法模拟研究   总被引：1，自引：0，他引：1       下载免费PDF全文

龙小刚  朱守彪 《地球物理学报》2015,58(7):2350-2365

台湾是菲律宾海板块与欧亚板块会聚、碰撞的产物,区内地质构造复杂,地震频发,变形强烈.为定量研究台湾地区变形特征及其动力学成因,文中运用二维不连续变形体弹性力学的有限单元计算方法,利用台湾地区1995—2005年GPS的观测结果作为边界约束,对台湾地区地壳变形场进行了模拟.计算结果显示,台湾中部地区主要呈压缩状态,但在台湾的东北部及南部地区出现了拉张的变形环境.变形程度最大的区域位于台湾东部海岸山脉及其附近海域,同时,计算给出台湾纵谷断裂滑移速率约为13.8~23.5mm/yr,由于纵谷断裂对变形的吸收,因此变形在纵谷断裂以西及西北地区迅速衰减.此外,计算结果还发现,计算给出的台湾岛上的速度值与GPS观测结果吻合得较好;计算给出的主应力方位与地应力观测及震源机制解结果也颇为一致.由此说明文中的有限元模型是合理的.此外,计算结果还表明,菲律宾海板块向北西方向的推挤,板块边界形状,观音、北康高地的抵阻,冲绳海槽扩张及琉球海沟的向南后撤以及断裂作用等,共同造成了台湾地区现今变形场的主要格局.  相似文献   

电导率各向异性的海洋电磁三维有限单元法正演   总被引：10，自引：8，他引：2       下载免费PDF全文

蔡红柱  熊彬  Michael Zhdanov 《地球物理学报》2015,58(8):2839-2850

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.  相似文献   

Multi‐method virtual electromagnetic experiments for developing suitable monitoring designs: A fictitious CO2 sequestration scenario in Northern Germany          下载免费PDF全文

Jana H. Börner  Feiyan Wang  Julia Weißflog  Matthias Bär  Ines Görz  Klaus Spitzer 《Geophysical Prospecting》2015,63(6):1430-1449

We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO₂ sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO₂ injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO₂ dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO₂, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO₂ plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.  相似文献   

Dynamic FE simulation of four‐story steel frame modeled by solid elements and its validation using results of full‐scale shake‐table test          下载免费PDF全文

Tomoshi Miyamura  Takuzo Yamashita  Hiroshi Akiba  Makoto Ohsaki 《地震工程与结构动力学》2015,44(9):1449-1469

Dynamic finite element analyses of a four‐story steel building frame modeled as a fine mesh of solid elements are performed using E‐Simulator, which is a parallel finite element analysis software package for precisely simulating collapse behaviors of civil and building structures. E‐Simulator is under development at the National Research Institute for Earth Science and Disaster Prevention (NIED), Japan. A full‐scale shake‐table test for a four‐story frame was conducted using E‐Defense at NIED, which is the largest shaking table in the world. A mesh of the entire structure of a four‐story frame with approximately 19 million degrees of freedom is constructed using solid elements. The density of the mesh is determined by referring to the results of elastic–plastic buckling analyses of a column of the frame using meshes of different densities. Therefore, the analysis model of the frame is well verified. Seismic response analyses under 60, 100, and 115% excitations of the JR Takatori record of the 1995 Hyogoken‐Nanbu earthquake are performed. Note that the simulation does not reproduce the collapse under the 100% excitation of the Takatori record in the E‐Defense test. Therefore, simulations for the 115% case are also performed. The results obtained by E‐Simulator are compared with those obtained by the E‐Defense full‐scale test in order to validate the results obtained by E‐Simulator. The shear forces and interstory drift angles of the first story obtained by the simulation and the test are in good agreement. Both the response of the entire frame and the local deformation as a result of elastic–plastic buckling are simulated simultaneously using E‐Simulator. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Finite‐element analysis of laminated rubber bearing of building frame under seismic excitation          下载免费PDF全文

Makoto Ohsaki  Tomoshi Miyamura  Masayuki Kohiyama  Takuzo Yamashita  Masashi Yamamoto  Naohiro Nakamura 《地震工程与结构动力学》2015,44(11):1881-1898

Finite element analysis is carried out for a building frame supported by laminated rubber bearings to simultaneously investigate global displacement and local stress responses under seismic excitation. The frame members and the rubber bearings are discretized into hexahedral solid elements with more than 3 million degrees of freedom. The material property of rubber is represented by the Ogden model, and the frame is assumed to remain in elastic range. It is shown that the time histories of non‐uniform stress distribution and rocking behavior of the rubber bearings under a frame subjected to seismic excitation can be successfully evaluated, and detailed responses of base and frame can be evaluated through large‐scale finite element analysis. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

应力扰动对地震周期和地震矩影响的数值模拟          下载免费PDF全文

翁辉辉  黄金水 《地震学报》2015,37(1):65-79

利用二维有限元数值模型,结合断层滑移弱化摩擦准则对断层滑动规律以及应力扰动对其影响进行了研究．数值计算结果表明,在均匀应力分布情况下, 平面断层滑动显示出典型的特征地震规律,断层面上的应力扰动对断层滑动规律产生影响,压应力增加明显延迟地震的发生时间,并增加地震释放的能量．应力扰动发生在地震破裂临界区时的影响比在震前滑移区时的影响显著．当发生在地震滑移区时,若应力扰动足够大,则压应力增大会造成地震发生时部分动力断层被暂时锁住,使得地震释放的能量变小,但可增加后续地震的能量; 而压应力减小则可导致地震规律产生更加复杂的变化,会即时触发地震．如果应力扰动发生在一个地震周期的早期,则触发的地震较小,但可导致随后的地震提前发生; 如果应力扰动发生在一个地震周期的后期,则会触发大地震．当应力扰动位于震前滑移区或破裂临界区时,小的扰动也可能产生类似的效果．应力扰动产生越晚,这种影响也越明显．应力扰动发生在破裂临界区的影响最明显．应力扰动的影响一般主要集中在应力发生扰动后的1—2个地震周期内．后续地震基本恢复无应力扰动时的特征地震规律．   相似文献   

three-dimensional sensitivity coefficients of apparent resistivity and preliminary application     

XIE Tao  LU Jun 《地震地质》2015,37(4):1125-1135

We calculate three-dimensional sensitivity coefficients distribution of apparent resistivity observation when Schlumberger array is used by using finite element method. Analysis results suggest that for the situation of one-dimensional positive or minus coefficient of surface medium, three-dimensional sensitivity coefficients distribution at surface shows similar patterns, and sensitivity coefficients distributions of different layered electric structures are also similar. There are two approximate ellipses at the two-dimensional surface plane between current electrodes and potential electrodes, where sensitivity coefficients are minus, and sensitivity coefficients at other areas are positive. Sensitivity coefficients at two approximate ellipses between current electrodes and potential electrodes are minus at the vertical section along monitoring line, while others are positive. From the three-dimensional view, minus sensitivity coefficients are at the two approximate half ellipsoids between current electrodes and potential electrodes when arrays are applied at surface. And coefficients near the electrodes are much greater than other areas. When resistivity of local areas at surface changes, we can qualitatively analyze the disturbing effects caused by the areas using three dimensional sensitivity coefficients distribution, and the analysis result can serve as reference for further experiment and numerical model quantitative analysis.  相似文献